SISSA 22/2009/EP 



O 
O 

o 



(N 



(N 



Comparison of two methods for the definition of the 

infinite time Hmit of 
the time evolution operator in scattering theory 



Gabor Zsolt Toth 



Q ' International School for Advanced Studies [SISSA ), 

cn ■ Via Beirut 2-4, 34OI4 Trieste, ItalyQ 



INFN, Sezione di Trieste, Italy ^ 



JlLi' Research Institute for Particle and Nuclear Physics, 



Hungarian Academy of Sciences, Pf. 49, 1525 Budapest, Hungary 

e-mail: tgzs@cs.elte.hu 



^ • Abstract 
^' 

Q>^ ' A comparison of the methods of averaging and adiabatic switching for the definition 

' of the infinite time hmit of the time evolution operator in scattering theory is pre- 

sented. The focus is on certain basic features which are relevant mainly in quantum 
^ , field theory. Specifically, the differences between the unitarity properties of the op- 

^ I erators obtained using the two methods, the difference between the necessary energy 
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1 Introduction 



The T — ^ ±cxD limits of the time evolution operator f/(0, — T) have an important role in 
scattering theory and also in some other areas of quantum physics. However, in general it 
is not completely obvious how to define these limits. In this paper we discuss certain basic 
aspects of this problem; we compare two methods for the definition of limT-^-too U{0, — T), 
which are mentioned in the literature as the method of averaging and the method of 
adiabatic switching. We emphasize certain features of the two methods and highlight 
certain significant differences between them. The discussion is done in the context of 
scattering theory, where the time evolution operator is used to produce the in and out 
states. 

In section [2] we introduce briefly the in and out states, the S-matrix elements and the 
time evolution operator. In sections [3] and H] we describe the method of averaging and the 
method of adiabatic switching, respectively. 

The main differences between the two methods are the following. 1) In general a 
renormalization of the free particle Hamiltonian operator is necessary, but the amounts 
of energy renormalizations that are needed are different for the two methods. 2) The 
limr^ioo ^(0, — T) operators obtained by averaging are in general not unitary, therefore 
normalization factors have to be included in the formulas which give the in and out states. 
Such normalization factors are not needed in the case of the adiabatic switching. 3) The 
Lippmann-Schwinger equations which hold in the case of the method of averaging do not 
hold in the case of the method of adiabatic switching. 

The difference between the two methods can be studied in perturbation theory as 
well. In particular, it can be seen that the two methods are identical to first order, but 
generally there are differences at higher orders. 

The differences that we discuss are relevant quite generally, in particular also for quan- 
tum field theory; nevertheless in quantum mechanical potential scattering the differences 
between the two methods usually vanish. 

In addition to the sections [21 [3] and H] described above, the paper contains section [5l 
where we present an illustration of the ideas described in sections [3] and H] in the case of 
the 0^ model, and section [6l where the conclusions are presented. Appendix A contains 
formulas relevant for section O 

We note that we were partially motivated to write this paper by recent development 
in scattering theory [I]-[l3], and that this paper is related to an earlier paper [U] of ours. 

2 The in and out states 

Let H be the total Hamiltonian operator that describes a certain scattering process. We 
assume that H takes the form H = Ho + fi'-f^int, where Hq is a Hamiltonian operator that 
describes free particles, gHi^^t is an interaction term and g is a coupling constant that can 
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be used as a variable for power series expansion. 

The S-matrix elements are scalar products of in and out states: 

Swv = {w , out\v , in) . (1) 

An important property of the in and out states is that they are eigenstates of H. 

In quantum mechanics standard formulas (see flbl [T6l [T7] ) that give the in and out 
states in terms of H, Hq and the eigenstates of Hq are the following: 

\v,in) = lim U{0,-T)\v), (2) 
\v,out) = lim U{0,T)\v). (3) 

T — ^oo 

In these formulas U{t2,ti) is the time evolution operator 

U{t2,ti) = e^^«*2e-*^(*^~*^)e"*^°*S (4) 

and \v) is an eigenstate of Hq. The in and out states \v,in) and \v,out) thus correspond 
to the eigenstate \v) of Hq, and limc,^o \v,in) = limc,^o \v,out) = \v). The energy of an 
eigenstate \v) will be denoted by E'^: Ho\v) = E^\v). The energy of \v,in) and \v,out) 
will be denoted by E^: H\v, in) = Ey\v, in), H\v, out) = Ey\v, out). The formulas ([2]) and 
imply that E^ = E^ 

In the following sections we aim to write ([2]) and ([3]) in a form that is valid generally, 
in particular also in quantum field theory. We shall see that certain modifications of ([2]) 
and ([3]) are necessary for this, and these modifications depend on the way the T ^ oo 
limit is defined. We will discuss two methods, which we mentioned in the introduction, 
for defining this limit. The equation will not be changed. 

We will write three different modified versions of ([2]) and ([3]); the first one is fl2Up and 
fl2T]) (or, with some regularizations explicitly denoted, fl2^ and fl5Ul) ). the second one is 
(148|) and (H9|) (or (l53i) and (15^ ). and the third one is (l55l) and (l56l) . The first version is 
based on the method of averaging, whereas the second and third versions are based on 
the method of adiabatic switching. We mention the third version only briefly, since it was 
the subject of [H]. 



3 The method of averaging 

In order to properly define the T — >■ oo limit of U{0, —T)\v) and U{0, T)\v) an averaging 
can be introduced in the following way (see e.g. [E]): 

UeiO,-T) = e [ dre-^^[/(0,-r), (5) 
Jo 

UMT) = e [ dre-^^f/(0,r), (6) 
Jo 
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and then 



f , in) 



lim lim ?7,(0, -T)|t;), 



(7) 



f , out) 



lim lim U,{{),T)\v). 
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The T — i> oo limit defined using this average is also known as the Abelian limit. 

The above definition for the in and out states is usually satisfactory in potential 
scattering (at least for sufficiently short range and regular potentials). However, in general 
(in quantum field theory, for example) one encounters divergences in the e +0 limit 
if one uses the definitions ([7j) and ([8]) and applies perturbation theory to calculate the 
components of the in and out states and the S-matrix elements. 

These divergences occur in any theory in which the energy of (the actual) |w,m) and 
|u, out) states, which are eigenstates of if, is shifted with respect to the energy of \v) by 
any nonzero amount. These shifts are called self-energy corrections. We use the notation 
5E^ for them (thus SE^ = — E^). We stress that the mentioned divergences are 
distinct from the ultraviolet divergences which often occur in quantum field theory, and 
also from the type of infrared divergences which are caused by soft or collinear massless 
particles in theories like quantum electrodynamics and quantum chromodynamics. We 
shall not discuss these types of divergences in this paper. In quantum field theory the 
divergences that we discuss are associated with disconnected vacuum bubble diagrams 
and with radiative corrections on external lines. A detailed explanation can be found e.g. 
in section 4.6 of [TB| . 

The problem of self-energy corrections can also be illustrated in a nonperturbative 
framework; for this purpose let us consider the case when the Hilbert space is finite 
dimensional and the spectrum of Hq is nondegenerate. Although the scattering is trivial 
in such a system, the self-energy corrections are generally nonzero. We consider the generic 
case when H is such that for sufficiently small but nonzero values of g all eigenvalues of 
H are shifted from the values that they take at g = 0. We assume that 1"^,^), for an 
eigenvector \v) of Hq, is defined by ([7j). At (? = we have \v,in) = \v), of course. For 
small values of g, the limits in ([7j) exist, but the shift of the eigenvalues of H has the 
consequence that {W\v,in) = for any eigenvector \W) of H, and thus \v,in) = 0. This 
can be seen in the following way. Let E^ be the eigenvalue for \W). Then we have 



which is an oscillating function of T since E^ ^ Ew- It is easy to verify that the averaging 
procedure yields zero for this function. This holds for all eigenvectors of H, and these 
eigenvectors form a basis in the Hilbert space, therefore the limits in ([7]) exist and |f , in) = 
0. This implies, in particular, that \v,in) is discontinuous at (7 = as a function of g. 



{W\U{0, -T)\v) = {W\v) exp[^(E° - Ew)T], 



(9) 
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3.1 Renormalization of Hq 



The standard solution to the problem described above is the renormalization of Hq, which 
consists in replacing the initial Hq operator in the definition of U(t2,ti) by a new one, 
while H remains unchanged. U{t2,ti) is thus replaced in ([5j) and (j6]) by 

f/-(t2,ti) = ei^rt2^iH^2-t^)^iHrti (10) 

and we have the formulas 

t/r"(0,-T) = e [ dTe-'^U''''{0,-T), (11) 
Jo 

?7f-(0,T) = e / dre-^"?7'^""(0,r), (12) 
Jo 



\v,in) = lim lim f/f'^(0,-T)|t;n, (13) 

e^+O T— >oo 

\v,out) = lim lim f/J^"^(0, T)|t;''^'^). (14) 

The vector \v'^'^^) appearing in fll3p and (IT^ should be an eigenvector of H^^^. We denote 
the eigenvalue for by H'q'"'\v'^'') = El''''\v''"^) . El'"' and AE^ = E'J''' - E° are 

called renormalized energy and energy renormalization, respectively. H^'^'^ has to have the 
same eigenvalue on |t>''™) as H has on \v,in) and \v,out), i.e. E^'^'^ = E^ has to hold. 

(More generally, one can replace the pair {Hq, H) by a suitable renormalized pair 
{Hq^^ , H^'^'^) , and in this case it is also possible that Hq = Hq^^. For the subject of this 
paper it has no significance whether H is changed or not.) 

In perturbation theory the condition that the e +0 limit in (fT3l) and (|T4|) should 
be convergent imposes restrictions on the Hq^'^ operator and on the \v'^^^) vectors at every 
order. 

In [IHl EOl EI] the following convenient choice is proposed: 

Hr = Ho + j dzAE, \v,){v,\ = J dzEr h){v,\ (15) 

and 

\vr) = h), (16) 

where the \vi) vectors are eigenvectors of Hq and constitute an orthonormal basis i.e. 

{vi\vj) = 6{i,j) , (17) 

where 6{i, j) is a Dirac-delta. The integration in f|T5l) is understood in a general sense; it in- 
cludes, for instance, summation over discrete states. The Dirac-delta is also a generalized 
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one. The index i labels all eigenvectors of Hq. The AEi are the energy renormalizations 
and E^'^'^ = + AEi are the renormalized energies. For all E^*^^ the equation 

Er = Ei (18) 

has to hold, where the Ei is the eigenvalue of the H operator on \vi, in) and \vi, out). The 
renormalized energies El^° can be calculated order by order in perturbation theory from 
the condition that the e +0 limit in (fT3|) and (|T^ should be convergent. 

The choice above for H^'^^ and {v^^'^) is not suitable in some cases (we refer the reader 
to section [23] for more details), nevertheless we conjecture that it is suitable for a large 
class of quantum field theories. 



3.2 Normalization of the in and out states 

The formulas ffT^ and ffn|) are still unsatisfactory because lim^^+o limT^oo ^^'^'^(0, — T) 
and lim^^+o liKiT^oo ^^^"^(0, T) are generally not unitary (more precisely, not scalar prod- 
uct preserving), even though U^'^'^{t2,ti) is unitary. This is caused by the property of the 
averaging procedure that it assigns in the T — ^ oo limit to functions which oscillate 
around 0. 

In order to illustrate the nonunitarity of lim^^+o lim-r^-i-oo ?7J'''^(0, — T) we consider 
again the case when the Hilbert space is finite dimensional and the spectrum of Hq is 
nondegenerate. We assume that H^'^^ takes the form fllSj) and g is sufficiently small so that 
the spectrum of H is also nondegenerate. Considering the definition of W'^^{0, —T) and 
the property that \v, in) is an eigenvector of H, we obviously have {v, in\U^'^'^{0, —T)\v) = 
{v,in\v). This implies that 

{v,in\v,in) = {v,in\ lim lim (/^^^{O, —T)\v) = {v,in\v). (19) 

However, {v,in\v,in) = {v,in\v) cannot hold if {v,in\v,in) = {v\v) unless \v,in) = \v), 
which is not the case generally. Therefore {v,in\v,in) ^ {v\v). In fact, the operator 
lim^^+o limT^oo ^^^"^(0, — T) projects \v) onto the (one dimensional) space spanned by 
\v,in), implying in particular that {v,in\v,in) < {v\v). 

On the basis of the above considerations, we can say that the correct form of ([7j) and 
dH]) should be 

\v,in) = lim lim ^ ^ U''''(0, -T)\v''''), (20) 



\v,out) = lim lim -,=i==f/r(0,T)|O, (21) 
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where 



/ renlfyrcn/Q -J^U^ron/Q _rp)\yTcn\ 

_ (t;-"|^f"(0,T)t[/f'^(0,T)|O 



The numbers Z^^^^m and Zy^^^out are real and are included in fl20l) and fl2Tl) in order to 
maintain the correct normalization of \v,in) and |f,OMt). 

We note that one could also consider Hq^^ operators which depend on e. Such a depen- 
dence can have an effect on the phase and normalization of lim^^+o hniT^oo U^'^'^iO, —T)\vi) 
and lime^+o limT^oo ^^"^(0, T)|f j). In particular, one can take e-dependent renormalized 
energies in (ITSj) . The requirement (|T8|) is relaxed then to 

hm = Ei. (24) 

It is not difficult to check that taking a renormalized energy level that depends on e as 

Er''{e) = Ei + iheci + 0{e^), (25) 

where Cj is a complex number, will change the normalization of 
lime^+o limT^oo ^^^"^(0, —T)\vi) by a factor 

^ (26) 



1 + Q 

and the normalization of lime^+o limy^oo U^^'^{0,T)\vi) by a factor 



(27) 



If Cj is real, then these factors are also real, but for a complex q they have a nonzero 
phase as well. In this paper we restrict ourselves to H^'^^ operators which are independent 
of e. 



3.3 Regularization and further remarks 

The expressions f l22l) and fl2^ are somewhat formal, since is usually not a normal- 

izable eigenvector, thus the denominator is not a finite number. If necessary, one can 
introduce a regularization to properly define these quotients. One possible regularization 
method is the following: one takes superpositions 

\vr,f^) = I da /,(^, «)!<-), (28) 
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where /i is a regularization parameter and /^(«, a) is a suitable function with the property 
that {vf'^,iJi\vY'^,ii) is finite for all i and j if /i 7^ 0. The original \vf'^) vectors should be 
recovered in the ^ limit. In the formulas (120|) . (12T]) . (l22l) and (l23l) the vectors \v''^^) 
should be replaced by 

Quantum field theories usually exhibit a type of divergence that is related to the 
infiniteness of the volume of the space. In particular, the powers of iJint are not well 
defined, because these powers involve vacuum bubble diagrams, and the product of the 
Dirac-deltas expressing momentum conservation in the vertices of a vacuum bubble di- 
agram is not well defined. In order to deal with this divergence one has to introduce a 
suitable regularization of the powers of the ifint operator, and thus of the time evolution 
operator. The parameter for this regularization will be denoted by A. i^Q™ has to be taken 
to depend on A, and it will not have a finite limit as A ^ 0, because the eigenvalues of 
H are divergent. All eigenvalues of H contain the same additive divergent part, therefore 
this type of divergence cancels out in the differences of the eigenvalues. An illustration 
will be given in section [51 

Taking into consideration the regularizations described above, the formulas (120|) . (1211) . 
and (125]) take the following form: 

\v,in) = lim hm hm hm ^==L==f/r(0,-T)|t;'--,/i),, (29) 
\v,out) = lim hin hm hm ^ [7^(0, T) (30) 

A^0 6^+0m-0T^oo ^Z„,A,e,/.,o«i(T) 



(^;^^^^ (0, -r)t;7f '^(O, -r) |^;'-°^ /i)^ 

(^;^^"^ (0, T)t[/f ^(0, T) |^;'^^^ /i)^ 

If one calculates a quantity involving in or out states, then generally the calculation should 
be done with nonzero A, e and /i, and the — > 0, e — > +0 and A — > limits should be 
taken at the end. For instance, one should calculate the matrix element {vi,in\0\vj,in) , 
where O is an operator, as limx^Q\im^^+Q\im^^Q{vi, in, e, fi\0\vj, in, e, fi) x. The regular- 
ization associated with A should be done on the whole product {vi,in,e, fi\0\vj,in,e, fi), 
not separately on \vj,in,e, jj) and \vi,in,e, jj) . This rule applies also when the S-matrix 
elements are calculated. 

It is worth noting that in quantum mechanical potential scattering (assuming short 
range and sufficiently regular potentials) the self-energy corrections are usually (thus 
the renormalization of Hq is not necessary) and the normalization factors Z^^in and Z^^out 
are usually equal to 1 after the T — > 00, n —>■ 0, e ^ +0 and A — >■ limits are taken 
(in fact, it is not necessary to introduce the regularization to which the A parameter 



Zv,\,e,^,in(^ T) 
Zv,\,e,^,out(T^ 
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belongs). We refer the reader to [H], section 3, example 2 for a hint about the type of 
argument that can be used to derive these results. In quantum field theory, however, the 
energy renormalizations and the normalization constants usually have nontrivial values. 
This can be seen in specific examples (e.g. in the 0^ model, see section [5]) by perturbative 
calculations. 

For a nondegenerate state |t>'''^'^) the phase of |w,m) defined by (|2U]) (or (12^ ) is such 
that {v'^^^\v^in) is real, and the same can be said of \v,out). This implies that \v,in) = 
1^,0^^) and thus {v,out\v,in) = {v'^^^\v'^^^) . In quantum field theory, in particular, the 
definitions fl29l) and fl30|) yield the result Sq-^ = 1 for the S-matrix element describing the 
vacuum-vacuum transition. The one-particle S-matrix element for a stable particle will 
be 5'ki;k2 = (ki|k2) = <5^(ki — ^2), where k denotes the momentum of the particle. 



3.4 Perturbation theory 

In this section we discuss some aspects of the application of perturbation theory to (l20l) 
and fl2Tl) . The perturbation series for \imT^oo{vi\Ue{0 , —T)\vj) and limT^oo{vi\Ue{0, Ty\vj) 
with unrenormalized Hq serve as basic ingredients for this. 

The perturbation series for \imT-,oo{vi\Ue{0, ~T)\vj) with unrenormalized Hq reads as 

lim {vi\U,{0, -T)\vj) = {vi\vj) + ^g- ^'^^ 



00 / .\ k 



^V^/ J P{ij) + eP{mk-ij)+e P{m2j) + eP{mij) + e 

(33) 



where the notation 



P{i3) = -{E^ - E°), (u) = {v,\^-^\v,) = {v.\H,^,\v,) (34) 



is used, and the integrations over the intermediate states are understood in a general sense, 
in particular they include summation over discrete states. Ef denotes the eigenvalues of 
Hq. The object (ij) is a generalized function of the variables i and j. 
For limT^oo(^i|^e(0, T)'''|f j) we have 

lta(.,|C/.(O.T)t|.,) = + ^9^^ 

00 / .\ fc 



J " ^ P{mit) + e Pim2t) + e P{mk-ii) + e P{ji) + e 

(35) 
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These formulas can be used as a basis for deriving perturbative formulas for the compo- 
nents of the in and out states, for the S-matrix elements and for other quantities. 

In order to take into account the renormalization of Hq one should replace in the 
formulas (l33i) and (135|) the P{ij) given in (IMI) by 



p'^-(zj) = ^(i?r-^r)' (36) 



and the (ij) given in by 



TT rrren 

{ijY^^ = {vl^'^'l ^lO- (37) 

9 

In this way one obtains series for \imT^oo{vl'^^\U^^'^{0, —T)\vj'^^) and 
\im'j'^oo{vr'^\U^'^^{0,Ty\vY'^) . It is assumed, of course, that the \vl'^^) vectors constitute a 
basis and their normalization is {vl'^'^lv^'^) = 6{i,j). Since Hq'^'^ depends on g, the series for 
limr^oo('yrl^/r''(0,-T)|vj'^") and limT-.oo{vr''\Ur''iO,Ty\vp) will not be power series 
in g; in order to obtain power series one has to expand their terms and collect the terms 
having the same power of g. 

We note that in the case when the Hilbert space is finite dimensional the Rayleigh- 
Schrodinger perturbation theory can be expected to be reproduced in effect. The integra- 
tions in fl55]) and fl5^ reduce to summations in this case. The divergences in the e —>■ +0 
limit without the renormalization arise from those terms in which a zero energy difference 
P{ab) occurs in a denominator. 

We consider now the case when an eigenvector \vi) of Hq belongs to a degenerate 
eigenspace. Let us assume that ifg*^" takes the form (fTSl) . We find at first order that \vi) 
must be such that (ji) = for all j ^ i for which Ej = E^. This condition is the same as 
the condition that one finds in Rayleigh-Schrodinger perturbation theory at first order, 
and it is always possible to choose the basis vectors in the eigenspace of Hq belonging to 
Ef in such a way that this condition be satisfied. At second order we find that 



E 



(jm) (mi) 
El - E? 



(38) 



has to hold for all values oi j ^ i for which E^- = E^, where the summation has to be 
done over those values of m for which E'^ ^ E^. This condition is stronger than the 
one in Rayleigh-Schrodinger perturbation theory; in the latter case (l38l) is required to 
hold only for those values of j for which \vj,in) remains degenerate with \vi,in) at first 
order. While the weaker condition of the Rayleigh-Schrodinger perturbation theory can 
always be satisfied by choosing suitable basis vectors, the stronger condition (1551) cannot 
always be satisfied, without violating the condition at first order, if the degeneracy of the 
eigenspace belonging to Ef is broken at first order. This shows that the form f|T5|) for 
ifg™ is not suitable for all cases. 
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Finally, it follows from the formulas fl55]) and fl5^ that the vectors 
lim^^^oo t^e(0, —T)\vj) and lim-r^oo t^e(0, satisfy the Lippmann-Schwinger equations 

hm f/,(0, -T)\v,) = \v,) + llg-^-^—^H,^,\im 11,(0, -T)\v,), (39) 

hm UMT)\v,) = \v,) + \ H,^, ]im UMT)\v,), (40) 

where e > 0. Similar equations can be written with H^'^^; these equations take the form 

lim U!'''(0, -T)\v^ = \vr) + ^ ^ [a - HT) lim f>r(0, -T)WD. 

(41) 

lnnJr{0,T)\v,) = {vf-) + _ ^ (g-go^) ji^m^ ^f" (0, T) ■ (42) 

4 The method of adiabatic switching 

Another method that can be used to define the T ^ oo limit of U{0, —T)\v) and 
U{0,T)\v), instead of the averaging procedure described in section [3l is the adiabatic 
switching of the interaction. In this case the in and out states are given by 

\v,in) = lim lim U,{0, -T)\v), (43) 

e— >+0 T^oo 

\v,out) = lim lim U,{0,T)\v), (44) 



where 

ft2 



Ueit2,ti) = Texp 

and 



(45) 



i^int,.(t) = e-^l*le*^«*/7i„te-^^°*. (46) 

In ( H5!) T denotes time ordering. In potential scattering this definition, like ([7]) and ([8]), 
is usually satisfactory. In general, however, lim^^oo f^e(0, — T) and lim^^oo f^e(0, T) are 
not convergent as e ^ +0, and, correspondingly, the vectors \v,in) and \v,out) cannot be 
defined. Similarly to the case of the method of averaging, this divergence can be related 
to the existence of nonzero self-energy corrections. 
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4.1 Renormalization of Hq 

The problem of divergence in the e +0 limit can be handled by the renormalization of 
Ho, i.e. by a suitable choice for H^'^"^, which is used to define U^^^{t2,ti) by substituting 

e-^\%i^rt(^H - HDe-^n'^^rt (47) 

for gHi^t,e(t) in (H5i) . The definitions for \v,in) and \v,out) are thus 

\v,in) = lim lim f/f"^(0, -T)|0, (48) 

|w,OMt) = lim lim f/f'^(0,T)|t;''^'^). (49) 

The vector |t>''™) appearing in fj48l) and fH9|) should be an eigenvector of ifg'^". The choice 
of i^p'^" is restricted by the requirement that the e — > +0 limits in these formulas should 
exist. 

In general, the Hq*^^ operators that are suitable for the method of adiabatic switching 
are different from those that are suitable for the method of averaging. In particular, 
it is not true that H^^^ has to have the same eigenvalues as H. In order to illustrate 
this point we consider the case when the Hilbert space is finite dimensional and Hq has 
nondegenerate spectrum. The Hq'^'^ operator is taken to be diagonal, i.e. of the form (|T5l) . 
It is not difficult to check, using the perturbation formulas fl58l) and (1591) below, that if 
the adiabatic switching is applied, then the necessary shift of the eigenvalues of Hq is 
not the full self-energy correction (i.e. the difference between the energy of \v,in) and the 
energy of |f )). Using the notation 6Ej = gSEj + g'^SEj + g^SEj + . . . for the self-energy 
correction (for a state \vj)) and AEj = g^Ej + g'^AEj + g^AEj + . . . for the required 
energy renormalization, we find that 

AEj = 5E] and AE| = (50) 

i.e. at first order the two quantities agree, but at second order the required energy renor- 
malization is only the half of the self-energy correction. For AE^ we obtain 

3^ (jm2)(m2mi)(mij) ^ {jm){mj){{mm) - (jj)) 

^ 3(E^J,^-K0)(E^^-E0)^^. 2(^^-^0)2 ' > 

where the notation of section l474l is used. The expression for 5E^, that is well known from 
Rayleigh-Schrodinger perturbation theory, is 

3^ (jm2)(m2mi)(mij) >^ {jm){mj){{mm) - (jj)) 

^ iK, - E'^KK, - E^) {E^-E'^y ■ ^^^^ 

This shows that there is no simple relation between AE^ and 5E^. 

Although the energy renormalizations are generally not equal to the self-energy cor- 
rections, the latter can of course be obtained from the eigenvalue equations for the in and 
out states. 
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4.2 Normalization of the in and out states, regularization 



Another significant difference from the method of averaging is that limr_+oo f/e™(0, — T), 
limy^oo ^^£^^"(0, T) and their limits as e — +0 are unitary (more precisely, scalar product 
preserving), therefore it is not necessary to include any normalization factor in fH3l) and 



If one takes renormalized energy levels in (1151) that depend on e, then the phase of the 
in and out states can change. In this paper we restrict ourselves to H^'^'^ operators that 
do not depend on e. 

Regularizations similar to those described in section 13.31 might be needed also in the 
case of the method of adiabatic switching. Taking into account these regularizations, the 
formulas (HSl) and (HH) have the form 



\v,in) = lim lim lim lim f/f'^(0,-T)|w'^°°,/i)A, (53) 
\v,out) = lim lim lim lim f/f'^(0,T)|^;■'^^/i)A. (54) 



4.3 Modification by phase factors 

In quantum field theory the use of the method of adiabatic switching seems to be restricted 
mainly to the vacuum state, and in this case the Gell-Mann-Low formula [221 [23] is used, 
which does not require the renormalization of Hq. 

In [13] we proposed the following formulas, which can be regarded as generalizations 
of the Gell-Mann-Low formula: 

\v,in) = lim lim n„,,(-T)[/,(0,-T)|t;), (55) 
\v,out) = lim lim n„,,(T)?7,(0,T)|w), (56) 

where li^^^iT) is a complex number of absolute value 1, i.e. a phase factor, given by 



n.,.(T) = 





\U,{TMv){v\ 


\Ue{^.T)\v) 


{v\ 


|f/,(0,r)|t;) 



(57) 



It may, of course, be necessary to introduce regularizations similar to those described in 
section [3l3l The phase factors in fl55|) and fl56|) have no effect on transition probabilities. In 
[H] we stated the conjecture that in a large class of quantum field theories the divergence 
in the e — > +0 limit is entirely canceled out by the phase factors, without any need for the 
renormalization of Hq. Nevertheless, the renormalization of Hq is needed in some cases 
also with (155|) and (156|) . 
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4.4 Perturbation theory 



The perturbation series for \imT^oo{vi\Ue{0, —T)\vj) with unrenormahzed Hq reads as 

_ _ ^ 

'v,\UjO, -T)\v,) = (v,\v,) + -^Q ^}'^J + V ( ^ 1 



^imA,Vi\Ue{,^, -T)\vj) = {vi\vj) + '^9 p^^-j-^ ^ ^ + /. \Jl ) ^'^ J '^"^i'^"^2 • • • dm^-i 

(58 



k=2 

{inik-i) {mk-imk-2) (mami) (mij) 



P{ij) + keP{mk-i]) + {k-l)e"' Pim^j) + 2tP{mi]) + e' 



where the notation fl34l) is used. This formula differs from fl33|l in the coefficients of e in 
the denominators. To first order, however, there is no difference between the method of 
adiabatic switching and the method of averaging. This imphes, in particular, that the Z 
normalization constants in fl2U]) and fl?I]) can be taken to be 1 to first order. 
For YimT^oo{vi\Uf_{T ,Q)\vj) we have 



Yim{vi\U^{T, {^)\vj) = {vi\vj) + -^9 p^j-^ _^ ^ + ^ \~h^ ) J ^"^i'^"^^ • • • dm^^i 

• (59) 



k=2 

(imi) {17111712) (mfc_2mfc-i) {rrik-ij) 

P{mii) + e P{m2i) + 2e"' P{mk-ii) + {k - l)e P{ji) + ke ' 



These formulas can be used as a basis for deriving perturbative formulas for the compo- 
nents of the in and out states, for the S-matrix elements and for other quantities. The 
renormalization of Hq should be taken into consideration as described in section 13.41 



Concerning the suitability of f[T5|l and the choice of the basis vectors in degenerate 
eigenspaces of Hq in the case of finite dimensional Hilbert spaces, one finds similar results 
as in the case of averaging. 

A remarkable consequence of the fact that the coefficients of e in the denominators of 
( 158|) . (159|) and in (l33l) . (135|) . respectively, are different, is that the vectors 
limT^oo Ue{0, —T)\vj) and limr^oo f^e(0, T)|t>j) do not satisfy Lippmann-Schwinger equa- 
tions similar to (l39l) and fHOl) . Similar statement can be made about 
limr^oo ^7^(0, -T) and limr^oo f/f ''(0, T) |ff . 



5 Illustration: 0^ model 

In this section we illustrate the ideas described in sections |3] and H] in the case of the 0^ 
model. We discuss the vacuum- vacuum S-matrix element Sq-q, the one-particle S-matrix 
element S'k^jka and the two-particle S-matrix element 5'k3k4;kik2- We calculate these S- 
matrix elements using perturbation theory to second order, and some third order terms 
will also be considered in the case of 5'k3k4;kik2- The results are 'bare', i.e. they are ex- 
pressed in terms of the initial unrenormahzed parameters of the full Hamiltonian operator 
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H . The (f)^ model is known to contain ultraviolet divergences. These are not relevant for 
the subject of this paper, nevertheless they can be avoided by introducing a momentum 
cutoff. The calculations are straightforward, but rather lengthy, therefore we shall not 
write many details of them. Basic formulas concerning the (j)'^ model, which are relevant 
for the calculations, are presented in Appendix A. We take H^^ to be of the form (ITSl) . 



5.1 The method of averaging 

5.1.1 The vacuum- vacuum S-matrix element 

Without energy renormalization we obtain to second order in g that 
lim (fi|?7e(0,T)t?7,(0,-T)|r])A = 
= 1+(Z1\\^ [ d'k^d^kgd^cd^kz, 2[6l(kA + k^ + kg + k^)]^ 

+ (—Yg^ f d'kAd^Bd^cd^D [6l(kA + kB + kc + kg)]^ 

\h ) J 4! UAi^B^Jc^^D [j-i^A + ujb + ^^c + ^d) + e]^ ' 

where r] = (27r)^2it~6- This quantity diverges as 1/e in the e — +0 limit and as in 
the A — >^ limit. (The A parameter is explained at the end of Appendix A.) 

As explained in section [31 energy renormalization and normalization factors are needed 
to turn (1601) into the correct vacuum- vacuum S-marix element. The energy renormaliza- 
tion for \VL) is denoted as Ai?n = gAE^ + g^AE^ + . . . . We restrict ourselves to a AEq 
that does not depend on e. It is straightforward to see that AE^^ has to be in order to get 
a finite result for limT^oo(fi|f/f °(0, T)tf/f =^(0, -T)\n)x in the e ^ +0 limit. Calculating 
to second order we obtain then that 

lim(^]|f/;-(o,T)t^7^(o,-^)|^]);,= lim(^]|f/,(o,T)tf/,(o,-^)|^]), + ^?2 — ^. (ei) 

T — >oo T— >oo fl t 

From comparison with flBUIl one can see that AE"^ has to be chosen as 

AEUX) = - [ ^'^Ad'^Bd'kcd^ko [5l{kA + k^ + kc + 

^ J 4! uaujb^c^d uja + u^b + + 

in order for limr^oo(^^|?/f ""(0, T)tf/f °(0, -T)|fi)A to be finite when e +0. AE^(A) 
diverges as 1/A^ in the A — > limit. 

For the constants 2'^ ^ a m and ^ a out "we obtain 



, . -i\ 2 I d^k^d^ksd^kcd^kz, r]' 

^Q,t,\,in — -L + 1 — ) 9 



h J J 4! UJAUJBi^C^D 
[SlikA + kB+kc + kp)]^ 

[j:{uJA + ujb + u!c + wd) + e][-j:{iOA + u}B + uJc + ujd) + e] 



(63) 
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and 2'n,e,A,m = Zfi^^ x,out- These quantities are finite as e ^ +0 but diverge as in the 
A ^ hmit. We note that the energy renormahzation has no effect on the Z^^^^x i^ and 
constants to second order in g. 
Using the results (l60l) - (l63l) in the formulas ([T]) and (l29l) - (l32|) we find that the final result 
for the vacuum- vacuum S-matrix element to second order in g is Sn;n = out\Q, in) = 1, 
in agreement with the remarks at the end of section I3.3[ 



5.1.2 The one-particle S-matrix element 

To second order in g the contributions to \imj'^oo{^2\Ue{0,TyU^{0, — T)|ki)A correspond 
to the graphs shown in figure [H The rate of divergence of the second order part is 1/e in 
the e +0 limit. 

In order to turn limT^oo(k2|?7e(0, T)^f/e(0, — T)|ki)A into the correct one-particle S- 
matrix element energy renormahzation and normalization factors are needed. From the 
requirement that the e +0 limit should be finite we obtain that the first order energy 
renormahzation AE^^ has to be zero. The same condition yields that the second order 
energy renormahzation has to be 



AEliX) 



d^k^d^kfi rj'^ 



3! UJAi^BOOcOo{uJA + UJb + i^c — ^) 

d^k^d^kfi r/2 



3! UaUJb^C^{^^A + + + ^) 

+AE'nW, (64) 

where k^ = k — k^ — k^ in Uc and uj = u(k). 

AEl{X) diverges as A — > 0, but the difference AEl{X) — AE^{X) has a finite limit (in 
fact, it is a constant function of A). 

For the constants ^k,e,A,m and .^k,e,A,oMt we obtain 



^k,e,A,m 



■^k,e,A,OMt — Zq ^ X + 



-I 



^ d^k^d^ks T]'^ 



,2 



h J J 3! uaujb^c^ 
1 



[jliuJA + ujb + ujc - + + + cuc - t^) + e] 

2 f d^k^d^kfi 1]"^ 



''^ \ h J ^ J 3! ujaoob^c^ 

[jliujA + ujb + ujc + + e][-j:{ujA + ujb + uJc + uj) + e]' 

This quantity is finite in the e +0 limit and the energy renormahzation does not 
have any effect on it. After these normalization factors and the energy renormahzation 
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are taken into account we find that the resuh for the one-particle S-matrix element is 
(k2, out\ki, in) = 6^{k2 — ki), in accordance with the remarks at the end of section 

The energy renormahzation lim\^o{^E'^{^) ~ ^Eq{X)) is equal to the energy correc- 
tion that can be obtained in the Lehmann-Symanzik-Zimmermann (LSZ) formalism. We 
note that both Zk^^^x,in/Zn,e,\,in Ai?^(A) — AE^{\) are independent of A. 



5.1.3 The two-particle S-matrix element 

In this section we consider the two-particle S-matrix element (ks, k4, 0Mt|ki, k2, m). The 
diagrams corresponding to the contributions to limT^oo(k3, k4|?7e(0, T)"'"?7e(0, — T)|ki,k2)A 
to second order are shown in figure [2j There are additional graphs to be included in the 
calculation which are obtained by interchanging and k^ in [2^, [2li, [2^, [2]F and [2^; these 
are not shown in order to save space, and we will not mention them explicitly in the 
following. Those terms which are divergent in the e +0 limit correspond to the graphs 
12^, [2][ and [2^. The rate of divergence is 1/e. 

We find that the energy renormahzation at second order is 

Ai?k^k.(A) = AEUX) + [AElW - AEUX)] + [AEl{X) - AE^(A)], (66) 

i.e. the energy renormahzation for the two-particle state is the sum of the energy renormal- 
izations for the constituent particles, if the vacuum energy renormahzation is subtracted. 
This is in agreement with the natural expectation. It is also natural to expect that one 
would find a similar result to all orders of perturbation theory and for any number of 
particles. 

For the normalization constants ^kik2,e,A,m and ^kik2,e,A,ojif we find that 

Zkik2,e,\,in = Zkik2,e,\,out (67) 

and ^kik2,e,A,m is given by the relation 

Zkik2,e, X,in Zki^e,\,in Zk2,e,\,in 



-'Q,e,X,in Zfi^^ X ifi Zq^^ x iji 



{6t 



We conjecture that this relation holds to all orders of perturbation theory, and a similar 
relation holds for states containing arbitrary number of particles. 

The joint effect of the energy renormahzation and the normalization constants is that 
the contributions corresponding to the disconnected graphs [2^, [2]F and [2^ are entirely 
canceled out from the expression for the two-particle S-matrix, and only the terms corre- 
sponding to the graphs [2^, [2]d, ^ and[2]i remain. 

At third order we restrict ourselves to the brief discussion of a certain type of graph, 
shown in figure[31 which describes interaction between the two particles and which contains 
a radiative correction on an external line. The terms corresponding to this graph are 
present in the expression for limT^oo(k3, k4|?7e(0, Tyu^{0, — T)|ki, k2)A, and some of these 
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a b c 



Figure 1: Graphs for the one-particle S-matrix element to second order 




e f g 

Figure 2: Graphs for the two-particle S-matrix elements to second order 

terms are divergent as 1/e in the e — > +0 limit. The energy renormahzation, in particular 
the modification of the matrix elements of i^int in (15^ and (15^ due to the second order 
part of the energy renormahzation, eliminates this divergence. 

At third order the modification of the denominators in fl33|l and fl35l) due to the 
energy renormahzation is also relevant. This has the effect that the Dirac-delta factor 
5{E^^ + E^^ — E^^ —E^^) in the contribution corresponding to the graph shown in figure [2]d 
that expresses energy conservation becomes 5(-Ek3 +E\^^ — -Eki ~E\^2)i i-^- the unperturbed 
energy levels are replaced by the proper ones. 
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Figure 3: A third order graph containing a radiative correction on an external line 



5.2 The method of adiabatic switching 

In order to get the vacuum-vacuum S-matrix element we start with the matrix element 
\imT^oo{^\Ue{0,TyU^{0, —T)\Q)x, for which we obtain to second order in g that 

lim (fi|f/,(0,T)tf/,(0,-T)|fi)A = 

^ ^ ^2 f d^k^d^ksd^kcd^k^ 7^2 [^3(k^ + kB + kc + ko)] 



h J J 4! uaujb^^c^d 4ii^A + ujb + ujc + ^d) + e] 

I 



2 



^ , _ I 2 f d^^k^d^kgd^^kcd^kz, r,^ [dlikA + kg + kg + kD)f 

This differs from f l60l) . obtained in the case of the method of averaging, in a factor of 2 in 
the first integral on the right hand side. 
Similarly to (l6Ti) . we have 

•OA p2 

hm {n\Ur{0, T)tf/--(0, -T)\n), = hm (fi|f/,(0, T)tf/,(0, -T)\n), + g\ ^. (70) 

The second term on the right hand side is of the same form as in fl6ip . As a consequence, 
the second order energy renormalization needed in the case of adiabatic switching is 1/2 of 
that needed for the case of the averaging, in accordance with the discussion in section I4.1[ 
Similar result can be found also for the one- and two-particle states. As we mentioned in 
section 14. 2[ normalization constants for the in and out states are not needed. The final 
results for the vacuum- vacuum, one-particle and two-particle S-matrix elements to second 
order are the same as in the case of the averaging. 

Similarly to the case of averaging, the contribution of the graph [3] to 
limT^oo(^|t^e(0, T)'''t/e(0, — T)|f2)A contains terms which are divergent as 1/e in the e — * 
+0 limit, but the energy renormalization, in particular the modification of the matrix 
elements of i^int due to the second order part of the energy renormalization, eliminates 
this divergence. 

As regards the formulas fl55|) and fl56l) for the in and out states, we refer the reader 
to [H] for their discussion and their application to the 0^ model. In the case of these 
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formulas those are the factors 11^, ^ (see flFTI) ) that have the same effect as the energy 
renormahzation has in the case of fHHl) and fH^ . 

6 Conclusion 

We discussed the methods of averaging and adiabatic switching for the definition of the 
infinite time hmit of the time evolution operator U{0, ±T) in scattering theory. We 
described three significant differences between the two methods. The first difference is that 
the operator obtained using the method of averaging is not unitary in general, therefore 
normalization factors are needed in the definition of the in and out states in order to 
ensure the conservation of probabilities. In the case of the method of adiabatic switching 
such a problem of unitarity does not arise. The second difference is that in the case 
of the method of averaging the energy renormalizations needed in the free Hamiltonian 
operator are equal to the energy shifts caused by the interaction, whereas in the case of the 
adiabatic switching these quantities are generally different. The third difference between 
the two methods is that Lippmann-Schwinger equations exist in the case of the method 
of averaging, but similar equations do not hold in the case of the method of adiabatic 
switching. At the level of perturbation series these differences are encoded entirely in 
the (seemingly minor) difference between the coefficients of e in the formulas (133|) . (l35l) 
and (158|) . (l59l) . We presented an illustration of the application of the formulas for the 
in and out states in the case of the 0^ model using perturbation theory to low orders. 
The final results that we obtained for the S-matrix elements were the same for the two 
methods, and they are in agreement with the expectations. It would be interesting to 
extend the illustration to all orders of perturbation theory and to other models. It would 
also be interesting to find general formulas for the energy renormalizations in the case of 
the method of adiabatic switching. 

Appendix A 

The 0^ model has the Hamiltonian operator H = Hq + gH[^t, where is the Hamiltonian 
operator of the free scalar boson of mass mo: 




(71) 



and the interaction term i^int is 




(72) 



The free boson field is 



0C 




d^k 



(73) 



(2v/?)3 



20 



where 



u;{k) = ^Jml + k^. (74) 
The particle annihilation operators a(k) satisfy the commutation relation 

Kk),a(k')t]=<5^(k-k')-^. (75) 
For the boson field we have the commutation relation 



X 



,t),dt(l>ix',t)]=t6\x-x'). (76) 



The multi-particle eigenstates of Hq are 

n 

|ki,k2,...,k„) = {l[^^a{k,y)\fl), (77) 



i=l 



where \Q) is the vacuum state of Hq. These states have the normahzation 

{Q\Q) = 1 (78) 

and 

(ki,k2, . . . ,k,|k;,k^, ...X) = - kp(i))'^'(k2 - kp(2)) . ..S%k^ - kp(„)), (79) 

p 

where the summation is done over the permutations of n elements. 
The eigenvalues oi Hq for |ki, k2, . . . , k„) and |Q) are 

^kik,...k„ - '^(ki) + a;(k2) + • • • + u;{K) (80) 

and 

E'n = 0. (81) 
It is straightforward to calculate the matrix elements of i?int; a few examples read as 

(ki,k2,k3,k4| / d^x :04:|n) = (27r)353(ki+k2+k3 + k4)-^— ^= (82) 

(ki,k2,k3| / d^x :04.|k4^ = (27r)353(ki+k2 + k3-k4)-^— ^= (83) 

(ki,k2,k3,k4,k5| y d^x : 0^ : Ike) = (ki, k2, k3, k4| y" d^x : 0^ : P>,b\k^-k^) 

+(ki, k2, k3, ksl y" d^x : 0^ : |Q)5^(k4 - kg) + ... , (84) 
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where uji = uj{ki). 

The regularization that is needed because of the divergence associated with the in- 
finiteness of the volume of the space and that is mentioned in section 13.31 is introduced 
in the following way: if a diagram contains vacuum bubbles, then the Dirac-deltas ex- 
pressing momentum conservation at the vertices in these vacuum bubbles are replaced by 
the regularized Dirac-deltas 5f(x) = Sx{xi)S\{x2)Sx{x3), where S\{x) = if x < —A/2 or 
X > A/2 and Sxix) = 1/A if -A/2 <x< A/2. 
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